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ABSTRACT 

We have carried out an extremely long integration-time (9000s beam" ^) 21-cm blind 
survey of 60 square degrees in Centaurus using the Parkes multibeam system. We find 
O \ that the noise continues to fall as ^Jtobs throughout, enabling us to reach an Hi column- 



^ , density limit of 4.2 x 10^^ cm"-^ for galaxies with a velocity width of 200 kms~^ in 

^ ' the central 32 square degree region, making this the deepest survey to date in terms 

. . \ of column density sensitivity. The Hi data are complemented by very deep optical 



observations from digital stacking of multi-exposure UK Schmidt Telescope -R-band 
films, which reach an isophotal level of 26.5 _Rmagarcsec~^ (~ 27.5 _B magarcsec"^). 
^ . , 173 Hi sources have been found, 96 of which have been uniquely identified with optical 

■ counterparts in the overlap area. There is not a single source without an optical 

counterpart. Although we have not measured the column-densities directly, we have 
inferred them from the optical sizes of their counterparts. All appear to have a column- 
density of Nhi = 10^°-^^='='^ -^^. This is at least an order of magnitude above our 
sensitivity limit, with a scatter only marginally larger than the errors on Nhi- This 
needs explaining. If confirmed it means that Hi surveys will only find low surface 
brightness (LSB) galaxies with high Mhi/Lb- Gas-rich LSB galaxies with lower Hi 
mass to light ratios do not exist. The paucity of low column-density galaxies also 
implies that no significant population will be missed by the all-sky Hi surveys being 
carried out at Parkes and Jodrell Bank. 
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1 INTRODUCTION 

Most of our knowledge of galaxy populations derives from 
samples collected in optical surveys, where there are strong 
selection effects. For instance most of the galaxies we 
know about are barely brighter than the terrestrial sky, 
while galaxies which have lower surface-brightnesses may be 



* Now at BAE Systems 



severely under-represented. Attempts to overcome optical 
selection effects through the use of better detectors, larger 
telescopes, longer exposures and sophisticated image pro- 
cessing have met with partial success. This remains a very 
difficult process, however, and the subsequent corrections 
which have to be made to small number statistics are both 
large and controversial (e.g. Impey & Bothun 1997, Disney 
1999). 

Low surface brightness (LSB) galaxies could be sig- 
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Table 1. Equivalent central surface-brightnesses (in _B-band) for 
different column densities and different Hi mass-to-light ratios 
(calculated using the derivation in Appendix It can be seen 
that if Mui/Lg remains constant at the value for optically- 
selected galaxies of ~ 0.3, it is necessary to go to very low 
column-densities to find LSB galaxies. However, if Mhi/Lb rises 
as the surface-brightness falls, it is possible to reach low surface- 
brightnesses without requiring low column-densities. 
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nificant in a number of contexts. They might add to or 
even dominate the luminosity-density and/or mass-density 
of galaxies as a whole (e.g. Fukugita, Hogan, & Peebles 
1998). Additionally, they could be responsible for much of 
the metal line absorption seen in QSO spectra (Churchill & 
Le Brun 1998), while their presence in large numbers could 
be important to theories of galaxy formation and evolution. 

Contemporary views as to the global significance of LSB 
galaxies are wide ranging, partly as a result of conflicting 
and usually indirect evidence, but mainly as a consequence 
of conflicting interpretation (Impey & Bothun 1997; Davies, 
Impey, & Phillipps 1999). 

BUnd searches in the 21-cm neutral hydrogen line have 
long been considered an alternative to optical surveys for 
finding LSB disc galaxies (Disney 1976). The narrowness 
of the line means that cosmic expansion will discriminate 
between extra-galactic and local hydrogen, so reducing the 
local background to the instrumental level beyond a redshift 
of a few hundred kms~^. Unfortunately, blind 21-cm surveys 
have until recently been severely limited by technical consid- 
erations, in particular by the tiny areas and small velocity 
ranges that could be covered to any depth in a practical 
time. 

Recent technical advances, including high sensitivity 
multibeam receivers and powerful correlators, have allowed 
much more ambitious blind surveys to be carried out. All- 
sky surveys of the southern hemisphere from Parkes (the Hi 
Parkes All Sky Survey, HIPASS) and of the northern hemi- 
sphere from Jodrell Bank (the Hi Jodrell All Sky Survey, 
HIJASS) are currently being completed. These have integra- 
tion times of 450s and 350s per beam respectively (Staveley- 
Smith et al. 1996; Lang et al. 2003), giving very similar 
sensitivities to sources smaller than the beam. To supple- 
ment these shallow surveys we have carried out the much 
deeper HIDEEP survey (9000s beam" i) of a small area of 
sky (4° X 8°) out to a velocity of 12,700 kms"i. 

The main motive for HIDEEP was to reach previously 
inaccessible surface-brightness levels. In general one might 
expect surface-brightness to be correlated with Hi column- 
density. Indeed, if one assumes that gas and starlight are 
distributed over proportionate areas it is easy to show that 
( Appendix IXt: 
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a relationship tabulated in Table Q To reach LSB galax- 
ies {/lo ^ 24 Bmag arcsec~2) which have 'normal' amounts 
of Hi (e.g. Mhi/Lb ~ 0.3) requires an Hi survey reach- 
ing down to Nhi ^ 3 x lO'^^ cm"^^ while to reach really low 
surface-brightness objects with /lo ^ 26.5 Smagarcsec"^ re- 
quires an Hi survey ten times more sensitive. Unfortunately, 
such low column-densities can be reached only by very long 
integrations, irrespective of dish size (as diffraction-limited 
beams decrease in angular area exactly in inverse propor- 
tion to dish area). In single-dish surveys, the detection of 
sources smaller than the beam will depend only on the mass 
of Hi (Mhi) per velocity channel in the beam, and their 
Hi column-densities {Nhi, measured in cm"^) will be irrel- 
evant. As we shall demonstrate (Sectionj^below) our survey 
is peak-flux limited so that, for detection: 
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where Speak is the peak flux and a is the noise per channel 
per beam. As Speak oc (M^j/d^AV) and cr oc l/{D'^t^^^), 
where D is the diameter of the dish, d is the distance to the 
source, and tint is the integration time, this can be re- written 
as; 
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where ki is a constant. Thus the maximum distance for 
source detection will be: 
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and the maximum volume in which such sources {Mhi, AV) 
will be detected is 
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(where Qt is beam size in sterads, Ni, is the total number 
of beams in the survey and T is the total duration of the 
survey) . 

Thus the volume surveyed (and hence number of detec- 
tions) per unit time is 
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i.e. short integration times per beam are favoured in order 
to find the most sources. 

However, short integration times imply that sources will 
only be detected nearby, and if they are too close they will 
over-fill the beam, reducing the amount of Hi within it. In 
these circumstances, Equation El must be adapted to 
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where Nhi is the column density in cm ^ and fc2 is a con- 
stant. As Qb = {^/Df 

Nhi k 
AV>Vt 

where fc is a constant. Equation|H|is independent of D, dish 
diameter, and is a mandatory requirement for detection be- 
cause a source which cannot be detected when it fills the 
beam certainly cannot be detected when it does not. In 



The HIDEEP Survey 3 




Figure 1. Detectability of galaxies with different column densi- 
ties in a survey. We consider the case of a galaxy with mass Mui 
and velocity width AV. Panel (a): the galaxy just fills the beam 
{^beam = ^gal) and has a peak flux S close to the survey limit 
(dashed line in the cartoon spectrum on the right hand side). 
This galaxy will have the lowest detectable column density. In 
panel (b) the galaxy is smaller, does not fill the beam and has a 
higher column density. The higher column density compensates 
for the beam dilution, and the galaxy will still be detected, as 
shown on the right hand side. Panel (c): as (a) except that the 
distance d' > d. The galaxy does not fill the beam, but the low 
column density yields a peak flux S lower than the survey limit. 
This galaxy will not be detected. 



other words, short integration-time surveys are only sensi- 
tive to high column-density sources, irrespective of dish size, 
a limitation which is seldom acknowledged (e.g. Zwaan et al. 
1997). (See Appendix IbI for a full derivation of Equation |2j 

We now show that the volume in which a galaxy can be 
detected only depends on its peak flux, and is independent 
of the actual column density, as long as the peak flux is 
higher than the survey limit. 

First let us consider the case of a galaxy with mass Mhi, 
and velocity width AT/, which just flUs the beam (case (a) 
in Fig. 1). We will also assume that it is at the peak flux 
limit S of the survey. As the latter limit only depends on 
the flux per velocity channel we find that for this galaxy 

Mhi oc 5'(AV)d^ oc Nhi rl^i « Nhi d^d^ (9) 
or 

S • A\/ oc NnilAO^aai (10) 

where in this case ISD,gai ~ Af2i,eam- As our survey is peak- 
flux limited and Ai}beam is constant, we find that, at fixed 
AV^, the limiting column density is directly proportional to 
the limiting peak flux. In other words, a galaxy with a col- 
umn density lower than the limiting column density will also 
have a peak flux lower than the limiting peak flux, and can 
thus never be detected. 

A more compact galaxy with the same {Mhi , AV) at 
the same distance must necessarily have a higher column 
density (case (b) in Fig. 1). This increase in column den- 
sity compensates for the beam dilution factor r^^j/r^g^^ = 
Af2gai/ Afijeam , aud the galaxy will still be detected, as its 
peak flux S is higher than the survey limit. 

In contrast, if we take the galaxy from our original ex- 
ample (which just fills the beam at distance d) and put it at 



Table 2. Volume sampled for galaxies of different masses and 
column-densities in Mpc^, for a constant velocity width of 200 
kms~i 
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a larger distance d' (case (c) in Fig. 1), the peak flux drops 
by a factor (d' /d)'^ and will not meet the survey limit. 

The detectability of a galaxy is thus independent of its 
column density, provided that its peak flux is higher than 
the survey limit. A galaxy filling the beam at the peak flux 
survey limit will have the lowest detectable column density. 
At a flxed [Mhi, AV), galaxies with higher column densities 
(necessarily) do not fill the beam, but will be detected over 
the same volume, as they will exhibit identical peak fluxes. 
Similar galaxies at larger distances will not be detected, in- 
dependent of column density, as they will drop below the 
peak flux limit. 

Table |21 demonstrates this for galaxies of different 
masses and column-densities. The 10^* cm~2 galaxy is be- 
low the column-density limit and will never be detected 
in HIDEEP. All the other galaxies are above the column- 
density limit and are therefore detected out to the distances 
set by their peak fluxes. The volume over which a galaxy 
above the column-density limit will be detected is deter- 
mined solely by the peak flux of the galaxy. 

In summary, 21-cm surveys have two constraints: (a) 
a peak-fiux limit given by Equation |3] in which dish size 
(D) is a distinct advantage, and (b) a surface-density limit 
where dish size is irrelevant, but in which integration time 
per beam is all-important. In a search for high surface- 
density Nhi objects, short integration-times per beam are 
favoured; in a search for low column-density (and there- 
fore low surface-brightness) objects, long integrations are 
mandatory. (An alternative way of looking at this is to 
note that larger dishes project the same system noise on 
to smaller areas of sky, because their diffraction-limited 
beams are smaller, and so have to contend with higher ap- 
parent sky-noise). Thus HIDEEP and HIPASS are comple- 
mentary: HIPASS with its relatively short integration-time 
per beam (450s) picks up large numbers of high surface- 
density sources but is insensitive to objects below 1.6 x 10^^ 
cm"2 (^y ^ 200 kms^i) while HIDEEP with its very long 
integration-time per beam (9000s) is the first blind survey 
(see Table|^ capable of reaching much lower surface-density 
Nhi, and hence surface-brightness, limits 

Before the advent of the multibeam system at Parkes 
it was neither profitable nor practical to make the very long 
integrations required to reach low column density. Thus the 
limits quoted in such surveys for the total amount of cosmic 
Hi referred only to high column-density clouds and galax- 
ies, though this was rarely acknowledged. Their integration 
times per beam were mostly far too short to detect the lower 
surface-density features. HIDEEP (see Appendix |bJ| should 
and indeed does reach 

Nhi > 2.1 x lO'^'AV (11) 

~ 4 X 10^** cm"'^ for a typical galaxy with a velocity width 
of 200 kms-i. 
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In principle, therefore, we should be able to reach galax- 
ies with lower surface-brightnesses than any detectable be- 
fore either in Hi or in the optical (see Table 1). 

The Hi survey is complemented by deep optical obser- 
vations, reaching an isophotal level of 26.5 magarcsec"^ in 
the 7?-band, covering three-quarters of the survey region. For 
our analysis, we assume a value of Hq = 75 kms~^ Mpc~^ 
throughout. 



2 PREVIOUS 21-CM BLIND SURVEYS 

Despite technical difficulties, blind surveys have been car- 
ried out, often using special techniques (see Tabled for de- 
tails). For instance Shostak (1977) re-examined the signals 
in the 'off-beams' of an NRAO 300-foot survey, where the 
'on-beams' were pointed at bright, optically selected galax- 
ies. Latterly, Zwaan et al. (1997) and Schneider, Spitzak, & 
Rosenberg (1998) have used the Arecibo radio telescope to 
carry out deep surveys. The results of such blind surveys 
have generally proved negative in the sense that very few 
previously-uncatalogued galaxies or intergalactic gas clouds 
(IGCs) were detected. 

However, as shown in Table|21 such surveys did not have 
sufficient sensitivity to low column density (Nhi) gas to ac- 
tually detect such objects. Quite apart from the theoretical 
arguments leading to Equation |H] the column-density sensi- 
tivity of any survey can be estimated retrospectively from 
its sensitivity to unresolved sources because (see Appendix 

iV;?r(Ay)=4.5xlO^°(^^^^j xAV (12) 

where F^^j' and Al/""' are the integrated flux (in Jykms^^) 
and the velocity width (in kms~^) of the source in the sur- 
vey with the lowest value of Fhi/AV and 59 is the source 
size in arcmin. When considering large LSB galaxies which 
could fill the beam, we set 59 equal to the beam FWHM 
to find the sensitivity to such systems. That the left hand 
side of Eauation ll2l is independent of dish size can be seen 
from the right hand side as both the minimum flux (top) 
and the beam size (bottom) are inversely proportional to 
the dish area. An examination of the various survey source- 
limits shows excellent agreement between the column den- 
sity as calculated according to Appendix |n| and separately 
from applying Eauation ll2l to the various source lists. For in- 
stance if we look at the deep drift-scan survey carried out by 
Zwaan et al. (1997), the Arecibo Hi Strip Survey (AHISS), 
and examine those sources which lie within its main beam, 
we can use these to calculate the column-density limit of 
this survey for a typical velocity width of 200 kms~^. We 
find that this limit is: 

NniiARISS) > 3.5 x lO" (13) 

i.e. 7Vj|7 (AHISS) ~ 10^^'^ cm'^ - almost precisely the 
minimum value found by the VLA follow-up observations. In 
other words, the AHISS survey was a shallow survey capable 
only of picking up sources down to a few times 10^^ cm~^ 
with a velocity-width of ~ 200 kms"\ The AHISS results 
do not, therefore, rule out the presence of a population of 
low column- density galaxies. It should be noted that the 



5(T column- density sensitivity for AHISS given in Table |3 
is substantially lower than found in Equation 1131 This is 
consistent with the finding of Schneider et al. (1998) that 
the limit for AHISS is well above 5a. 

The other recent deep Arecibo survey, the Arecibo Slice 
(Schneider et al. 1998; Spitzak & Schneider 1998) operated 
rather differently. Patches of sky about 2 beam diameters 
apart were followed for 60s, and any sources tentatively 
picked up were then re-scanned using a grid search. That 
meant that the sensitivity varied by a factor of ~ 4 over 
the total 55 square degrees searched. The median flux of the 
sources detected over the velocity range 100 to 8340 kms~^ 
was 2.69 Jy kms^^ as against the median value for HIDEEP 
of 1.96 Jykms^^. Because of the short 60s integration time, 
the sensitivity to low-column-density sources was poorer 
than AHISS. Nevertheless there were interesting results. 
Only half the galaxies detected were in any optical cata- 
logue and about a third were LSB galaxies. The galaxies 
have much lower bulge-to-disc ratios than found in optically- 
selected samples. The median Mhi/Lb of the survey was 
O.89M©/I/0 and the galaxies with larger Mhi/Lb ratios 
had lower surface-brightnesses. 

Henning's (1992; 1995) survey with the Green Bank 
300-foot telescope was largely behind the galactic plane and 
therefore not comparable with HIDEEP. The Arecibo Dual 
Beam Survey (ADBS; Rosenberg & Schneider 2000; 2002) 
covered a larger area of sky than the other surveys but to 
a much shallower depth, and so could not set any limits to 
the population of low-column-density galaxies. 

Although blind 21-cm surveys are, in principle, the ideal 
way of circumventing optical selection effects and looking for 
LSB galaxies and IGCs, the weakness of the 21-cm signal 
and the system noise make it very difficult to find sources 
unless they have high column-densities. Only very long in- 
tegrations are capable of reaching low column-density lim- 
its, and possibly finding objects of lower surface-brightness 
than can be seen optically - thus the interest of HIDEEP. 
We claim that the arguments summarised in EquationQand 
TableQdemonstrate that we should, for the first time at 21- 
cm, be capable of locating such objects. Even so, the limita- 
tions of all such blind Hi surveys should constantly be kept 
in mind. Such surveys have lower sensitivity to broader-line 
sources (see below) and may, particularly at low column- 
densities (< 10^^ cm~'^), be severely affected by ionisation 
and spin temperature effects (Section |4j. Hi surveys, there- 
fore, can set only lower limits to the number of LSB galaxies 
and IGCs in the cosmos. 



3 THE HIDEEP SURVEY 
3.1 The Hi data 

HIDEEP was carried out in a region of Centaurus centred 
on a = 13''40'"00", 5 = -30°00'00" (J2000) with 1024 spec- 
tral channels covering -1280 to 12,700 kms"^ Due to the 
shape of the Parkes multibeam footprint, the survey has a 
sensitivity as good as or better than HIPASS over 6 by 10 
degrees with a uniform sensitivity over the central 4 by 8 
degree area. This region lies in the supergalactic plane, 30 
degrees from the galactic plane. The HIDEEP volume in- 
cludes the Gen A group (Banks et al. 1999) and the outer 
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Table 3. Blind Hi sureys 
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parts of the Centaurus cluster. The observations were car- 
ried out in the southern autumns of 1997, 1998, 1999, and 
2000. 

The data were processed using the standard multibeam 
reduction techniques, as described in detail by Barnes et 
al. (2001). Continuum sources were removed using luther 
(Wright & Stewart 2003, in preparation). Once integrated, 
the data take the form of a cube with voxels (3D pixels) 4 by 
4 arcminutes on a side and 13.2 kms~^ deep. The half-power 
beam width is 15 arcminutes and the data were smoothed 
in the velocity direction (to reduce ringing) as part of the 
reduction process, giving a velocity resolution of 18 kms~^. 
The data in adjacent voxels are therefore not entirely inde- 
pendent. 

The sky was Nyquist sampled 50 times and some 1800 
separate samples contributed to the signal in each voxel. 
Median filtering of this large sample greatly reduces inter- 
ference while the data were all taken at night to avoid solar 
radiation entering the beam sidelobes - a major source of 
noise during the day. 

The final data cube can be examined in 3 planes: {a, 5), 
{S,V) and {V,a), where V is the velocity direction. All 3 
planes are used for finding and measuring sources. Figure |21 
shows a {S, V) slice of the cube, showing the strong galac- 
tic signal at kms"^ as well as 12 other sources including 
NGC 5236 (M 83) at 500 kms'^ and GSC 7265 02190 at 
11875 kms"^ (Willmer et al. 1999). Continuum sources have 
been removed and the nature of the remaining noise, against 
which sources must be found, can be seen. The ripple seen 
just below -34° decl. is the residual of the strong continuum 
from the southern radio-lobe of IC 4296. The increased noise 
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Figure 2. A Decl. - Velocity slice through the HIDEEP cube. 
At least twelve sources can be seen, including Messier 83 at 500 
kms-i and GSC 7265 02190 at 11875 kms'^ (Willmer et al. 
1999) 

at the edges of the cube is the result of poorer sampling in 
these regions. 

Figure 13 demonstrates the effectiveness of long integra- 
tions. We plot the median noise in mjy beam" ^ channel"^ 
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Figure 3. Fall off of noise with integration time. Points are for 
HIPASS, the intermediate HIDEEP field with 5,600 seconds in- 
tegration time, and the full HIDEEP field with 9,000 seconds 
integration time. The line shows the theoretical 1 / \/ti^ fall-off 
normalised to the noise in HIPASS. 
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Figure 4. Fall off of noise with smoothing. The curves show 
^0.4 g^j^jj j-j^g theoretical curve from Poisson statistics, Af"'^ 
(where A'^ is the number of channels smoothed over). 



for integrations of 450 sbeam"^ (HIPASS), 5,600 sbeam"^ 
(Minchin 2001) and 9,000 sbeam"^ (HIDEEP). As can be 
seen it falls to 3.2 ±0.3 mjybeam"^ channel^^ against time 
in accordance with the theoretical 1 / yjtobs- 

Figure |1] shows that smoothing in the velocity plane is 
a much less effective way of reducing the noise. We have 



smoothed the data with a Hanning filter and removed every 
other channel from the smoothed cube in order to leave only 
the independent channels. This has been repeated to form 
cubes smoothed over 2, 4, 8, 16 and 32 channels. For white 
noise, the noise should fall as 1 / ^/ N^han (where N^han is the 
number of channels smoothed over), however it can be seen 
that this is not the case for Nchan > 4. Beyond this the fall 
off is shallower than predicted by Poisson statistics - closer 
to N'^j^^f^, although it is not well described by a power-law. 

3.2 Source finding 

The HIDEEP cube was searched three times in an inde- 
pendent manner. Two different people searched through the 
cube by eye, inspecting every channel and noting down the 
sources found, and the third search was carried out by an 
automated finding routine based on peak-flux detection and 
template fitting. This routine identified points higher than 
4.5-CT on a Hanning-smoothed cube and fitted Gaussian tem- 
plates with a large range of widths at these points, de- 
manding a correlation of better than 0.75. The three lists 
given by these searches were then compared and any sources 
found two or more times were accepted. The remaining dis- 
puted sources were examined by a third member of the team 
for a final decision. This team member had not previously 
searched the data cube. 

This gave us a final list of 173 sources, all of which 
have been judged as real by at least two members of the 
team acting independently. Accurate positions for all the 
sources were found by forming zeroth-moment maps around 
their positions and velocities and fitting Gaussians to these 
maps. This gives a positional accuracy, as judged from those 
sources that can be securely identified with optical galaxies, 
of around 2' (see Figure|5J. The Hi parameters of the galax- 
ies were measured using the mbspect routine in miriad 
which provides measures of the velocity width, the noise 
in the spectrum and the peak flux of the source as well as 
robust estimates of the integrated flux and systemic velocity. 

3.3 Completeness of HIDEEP 

The form of the selection present in the HIDEEP survey has 
been analysed by plotting the integrated flux of the galaxies 
against their velocity width (Figure |SJ. For a survey limited 
solely by the total flux of the galaxies, the selection limit 
on this graph would be a horizontal line. If the best-possible 
real-world selection was made, i.e. selection purely by signal- 
to-noise ratio using optimal smoothing in the velocity plane, 
then the selection limit would be a line with a slope of 0.5 
on a log- log plot (assuming SNR oc l/VAV", which is not 
strictly true as shown in Figure 2J, shown here as a dashed 
line. The solid line on the graph shows a selection limit based 
on peak flux {Fhi oc AV20); only 16 of the 173 sources fall 
below this line and only 8 of these are below by more than 
l-(j. Sources appear to fall further below the line at the low- 
flux, low-velocity-width end, but this is probably due to the 
larger errors in this region. The peak-flux selection limit is 
shown in Figure 1^] It can be seen that this explains well the 
selection limit of the survey. 

It can also be seen in both graphs that there is a paucity 
of galaxies narrower than approximately 4 channel-widths 
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Figure 5. Selection limits in velocity-width - integrated flux 
space. The theoretical Scr limit for selection based on constant 
signal-to-noise using optimal smoothing is shown by the dashed 
line and the 3(t limit for peak-flux selection {Ff[j/AV20 = con- 
stant) is shown by the solid line. 



Figure 6. Selection limits in velocity-width — peak-flux space. 
The 3(T (9.6 mjy) limit used above can be seen to be a good 
match to the selection limit of the data. The dashed line shows 
AV20 = 4 channels (52.8 kms~l). 



(52.8 kms^^). This further selection effect, thought to be 
due to galaxies narrower than this being indistinguishable 
from interference, is investigated in detail by Lang et al. 
(2003). 

The completeness of the HIDEEP catalogue has been 
calculated by looking at how the source counts vary with 
peak flux, as shown in Figures and |H1 For this analysis, 
41 sources that may have been detectable beyond the upper 
velocity limit were excluded in order to make a purely flux- 
limited sample. It can be seen from both flgures that the 
peak flux completeness limit is 18 mJy, or around 5.5a. This 
is around the level where the completeness limit is expected 
to fall. It can also be seen that large numbers of sources are 
found between the completeness limit at 18 mJy and the 
selection limit of ~ 10 mJy. 

3.4 The optical data 

To identify LSB galaxies it was necessary to obtain very deep 
optical data to compare with the radio. Accordingly eight 
1-hour 6° X 6° 7?-band Tech-Pan films were exposed at the 
UK Schmidt Telescope, centred on 13''39'"50" ,-30°00'12" 
(J2000). These were digitised, linearised and stacked using 
the SuperCOSMOS machine at the UK Astronomical Tech- 
nology Centre in Edinburgh (Hambly et al. 2001). The final 
image was then calibrated using the magnitudes of unsatu- 
rated ESO-LV (Lauberts & Valentijn, 1989) galaxies within 
the region, which yields a calibration accuracy of approxi- 
mately 0.2 magnitudes. The Tech-Pan films used go 1 mag 
deeper than the IIIaF plates previously used at the UKST 
(Parker & Malin 1999) and the digital stacking gives a fur- 
ther gain of over a magnitude compared to a single expo- 
sure (Schartzenberg, Phillipps, & Parker, 1996). The limit- 
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Figure 7. Completeness of the HIDEEP survey: Source count 
against peak flux (5pet,fe)- The histogram shows numbers found 

— 5/2 

in each bin of peak flux, the curve represents N{Speak) <^ ^peak 
as expected for a flux-limited survey. 

ing surface-brightness reached for small objects within the 
image is then 26.5 _Rmagarcsec~^, equivalent to between 27 
and 28 _Bmagarcsec~^. 

Within the overlap region between the radio and opti- 
cal images we find 96 Hi sources which have been uniquely 
identified with optical counterparts. 
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Figure 8. Completeness of the HIDEEP survey: Log-log plot of 
source count against peak flux. The solid line has a slope of —3/2, 
as expected for a flux-limited survey. 

Sources have been identified with galaxies on the op- 
tical image, firstly on the basis of positional coincidences. 
Multibeam survey positions are generally accurate to 2 ar- 
cminutes (Koribalski et al. 2003) , but this can be checked for 
the 65 optical counterparts which have previously published 
optical velocities or for which we have our own optical spec- 
troscopy or 21-cm interferometry data. A comparison can 
be made between the (radio - optical) offsets of these firm 
detections with the remainder (Figure lOJ. A Kolmogorov- 
Smirnov test confirms that there is no significant difference 
between the distributions, implying that the purely posi- 
tional coincidences can generally be trusted. There may still 
be one or two incorrect identifications, but the tail of off- 
sets out to 6.5 arcminutes probably reflects the positional 
accuracy of the HIDEEP survey. 

In the overlap area 59% of the sources are identifled 
with previously catalogued galaxies with matching redshifts, 
24% with previously catalogued galaxies without redshifts, 
and 17% are previously uncatalogued galaxies. It appears 
that there are no intergalactic gas clouds unassociated with 
optical counterparts: all those galaxies which have not been 
uniquely associated with a counterpart have more than one 
plausible candidate. 

We have measured effective radii and surface- 
brightnesses by fitting to the surface-brightness profiles of 
the optical sources using the ellipse routine in IRAF. The 
total magnitudes used were determined by SEXTRACTOR 
(Bertins & Arnouts 1996). The surface-brightness distribu- 
tion fFigure lTUl is much broader than one finds in optically 
selected samples. A Kolmogorov-Smirnov test of our distri- 
bution and the surface-brightness distribution of the ESO- 
LV shows that the hypothesis that both are drawn from 
the same parent population has a significance of less than 
1% - this is due to the larger number of LSB galaxies seen 
in the HIDEEP sample. This confirms that Hi surveys do. 




Offset (arcmin) 

Figure 9. Comparison between the cumulative distributions of 
offsets for the firm identifications (solid line) and the less certain 
counterparts (dashed line). 




i^.„ (R) 

Figure 10. Number of galaxies found in each surface-brightness 
bin. The solid line shows the distribution of observed surface- 
brightnesses, the dashed line shows the distribution of surface- 
brightnesses after correction for galactic absorption, cosmological 
dimming, and inclination. 



as expected, avoid some of the surface-brightness selection 
effects present in optical surveys. A fuller description and 
analysis of the optical properties will be given in a second 
paper (Minchin et al. 2003, in preparation). 
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Figure 11. Hi mass - velocity width relationsliip for 
HIDEEP sources. Tlie best fit to the HIDEEP data (A^ = 
0.42tJ5;57M^J*2='=° °^^) is shown by a solid line, while the fit of 

1/3 

AV = 0.15M^j (Briggs & Rao 1993) is shown as a dashed line. 
3.5 Internal Hi correlations 

A correlation between Hi mass and observed velocity width 
of the form AV oc M^j (where AV is the velocity width 
uncorrected for inclination) is expected in the HIDEEP data 
as it has been seen in optically-selected samples (e.g. Briggs 
& Rao 1993) and as it would be the consequence of the Hi 
TuUy-Fisher relationship. 

That such a relationship can be seen in the HIDEEP 
data is shown in Figure ITTl The solid lines shown here is for 
the best fit of AV20 = 0.42tE;;5°Af^f ^^^^^^ . The dashed 
line shows AV = 0.15Mj/j as found by Briggs & Rao (1993). 
The best-fit slope found for the HIDEEP data is 2a shal- 
lower than this, however this may be due to the selection 
against narrow velocity- width galaxies seen earlier. We can- 
not, therefore, conclusively say that our sample shows a dif- 
ferent relationship to that found by Briggs & Rao. 

It is also expected that there will be a relationship be- 
tween peak flux and the value of Fhi /AV20 - the peak flux 
of a top-hat function with a width of the 20% width of the 
source and containing the same total flux. This relationship 
is important for calculated Hi mass limits in a peak-flux lim- 
ited survey, as it is necessary for relating the peak-flux limit 
to the integrated flux (Fhi) on which the Hi mass depends. 

However, as this relationship depends on the profile 
shape it may well vary with Hi mass. This is investigated 
in Figure fT^ This figure shows that there is no dependence 
of the ratio J'^^j^^y Mhi- the slope of a fit to the points 
is statistically indistinguishable from zero. We therefore take 
the ratio to be a single number for all Hi masses, using the 
median value: 1.7±0.3. This value was used to calculate the 
3(T peak- flux limit in Figure |S] 

The Hi mass distribution of the HIDEEP sources is 
shown in Figure 1131 This shows both the distribution of 
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Figure 12. Comparison of the ratio p^'j^^v ™^th Mui- The 
slope of the best-fit to this data is indistinguishable from zero, 
we therefore use the median value, 1.7 it 0.3 to describe the ratio 
across all Hi masses. 

all the sources in HIDEEP (solid line), the distribution if 
the Centaurus A group galaxies are ignored (dashed line) 
and the distribution of galaxies included in the optical sam- 
ple (dotted line). It can be seen that most of the low-mass 
galaxies are, unsurprisingly, in the Centaurus A group but 
we find only one high-mass galaxy (Messier 83). There are 22 
hydrogen giants with Mhi > 10^°Mq, but only NGC 5291 
has a mass greater than 3 x 10^" M0. There are no galax- 
ies detected outside the range lO'^M© < Mhi < lO" Mq. 
The optical sample is similar in shape to the full sample, 
but contains fewer sources. This is mainly because it covers 
a smaller area. Some sources have also been omitted as a 
unique optical counterpart could not be identified. 

3.6 Sensitivity and comparisons with previous 
surveys 

Using the completeness limits and relationships found in 
the Hi data, we can calculate the sensitivity limit of the 
HIDEEP survey to galaxies of different masses and different 
column-densities. The mass (in solar masses) is related to 
the integrated flux by the equation 

Mhi = 2.356 x lO^fj/jdl/pc (14) 

(for Fhi in Jykms"^) and the total flux is related to the 
peak flux and the velocity width (Figure by 

Speaft = (1.7±0.3) X (15) 
AV20 

This allows the mass to be related to the peak flux as: 
Mhi ^ 1.386 x lo'^ V^feAVaodMpc (16) 

However the velocity width is a function of the mass, 
AV ~ 0.15M^''j'* (although with a large scatter). If we use 
the relationship found by Briggs & Rao (1993) then 
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Figure 13. Distribution of Hi masses. Tlie solid line shows the 
distribution for the whole survey, the dashed line shows the dis- 
tribution excluding galaxies in the Centaurus A group and the 
dotted line shows the distribution of galaxies included in the op- 
tical sample. 



Mhi ~ (0.15 X 1.386 X 10^ Vafc^Mpc)^^^ (17) 
which gives: 

Mhi ^ 3.0 x 10''5^/f,4fpc (18) 

and putting in the completeness limit of Speak ~ 0.018 Jy 



7.24 X 10^ 



Mpc 



(19) 



which is the sensitivity limit for sources of different masses at 
different distances. Since the relationship between velocity- 
width and mass has a large scatter this is not the absolute 
detection limit: sources narrower than predicted by this re- 
lationship will be seen further out, and sources wider than 
predicted will be found over a smaller volume. However, it 
does indicate the distance to which most sources of a given 
mass will be seen. 

The relationship between mass and velocity-width sug- 
gests the possibility of a selection effect in Hi surveys that 
appears to have been neglected in many previous surveys: 
that of a minimum believable velocity width. Below a cer- 
tain mass, the velocity-widths of many of the galaxies will be 
smaller than the minimum believable velocity-width of the 
survey - the width at which a peak in the data is recognised 
to be a galaxy rather than a noise peak - and will not be 
catalogued. A fuller analysis of this effect, using data from 
the Hi Jodrell All Sky Survey, is given in Lang et al. (2003) . 

Examination of the HIDEEP data shows the minimum 
believable velocity-width to be around 4 channels wide for 
AV20, or 52.8 kms~^. This selection effect means that de- 
tecting large numbers of low-mass galaxies will require not 
only sensitivity but also narrow channel- widths. However, 
this effect could remove at most ~ 40 per cent of 10^ A/0 
galaxies here. If the thermal broadening of the Hi is also 



Figure 14. Hi mass limits for three surveys: HIDEEP (solid line), 
the HIPASS Bright Galaxy Catalogue (BGC; dotted line) and the 
Arecibo Hi Strip Survey (AHISS; dashed line). 



taken into account, this percentage will fall. It is there- 
fore unlikely that this will significantly change the shape 
of Hi mass functions down to their current mass limits. 
For HIDEEP, the distance at which galaxies that would be 
above the mass limit will fall below the minimum believable 
velocity-width limit is about 18 Mpc. 

Figure 1141 compares the Hi mass sensitivity of three 
surveys: HIDEEP, the HIPASS Bright Galaxy Catalogue 
(Koribalski 2003; BGC) and the Arecibo Hi Strip Survey 
(AHISS). A minimum believable velocity- width of 4 chan- 
nels has been applied to all these surveys although galaxies 
narrower than this can be seen if they have high peak signal 
to noise ratios. This will apply particularly to the HIPASS 
BGC, where the peak-flux cutoff of 116 mjy is more than 
8cr. For AHISS a limit of 5.25 mJy, or seven times the stated 
noise value of 0.75 mJy per channel, was adopted. This ap- 
pears consistent with an analysis of the fiuxes and velocity 
widths from Zwaan (2000) and with the analysis of Schnei- 
der et al. (1998) which concluded that the AHISS survey was 
limited at 7cr due to the method of confirming sources. Only 
the primary beam area was used in calculating the volume 
covered by AHISS. 

Having estimated the sensitivity to sources in the 
mass-limited regime, we can also estimate the sensitivity 
to sources in the column-density limited regime. Column- 
density sensitivity has often been presented, in a similar 
manner to optical surface-brightness, as not having a de- 
pendence on distance. However the relationship between /S.V 
and Mhi and the dependence of column-density sensitivity 
of Al/ (see Eauation ll2t means that this is not the case for 
single-dish surveys - the tendency of higher mass sources 
to have larger velocity widths means that they will (a) be 
seen to greater distances due to their higher mass and (b) 
will have a higher column-density limit due to their larger 
velocity width. 
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Figure 15. Hi column-density limits for three surveys: HIDEEP 
(solid line), the HIPASS Bright Galaxy Catalogue (BGC; dotted 
line) and the Arecibo Hi Strip Survey (AHISS; dashed line). It 
can be seen that HIDEEP is considerably more sensitive than 
previous surveys to low column-density sources. 



Figure 16. Coverage of Muj - Njji space for three surveys: 
HIDEEP (solid line), the HIPASS Bright Galaxy Catalogue 
(BGC; dotted line) and the Arecibo Hi Strip Survey (AHISS; 
dashed line). This shows that neither the HIPASS BGC nor the 
AHISS were likely to find low surface-brightness giant galaxies. 
However, if these exist and are not very rare, they should be found 
by HIDEEP. 



The column-density of a source filling the beam (in 
Mq pc~^) is given by: 

where 6 is the HWHM of the telescope beam in radians and 
dpc is the distance of the source in parsecs. It can easily be 
seen that Nhi is only independent of distance if the sensi- 
tivity of the telescope to Mhi goes as d? . This would only be 
the case if galaxies at all different Hi masses had the same 
velocity-width or if the survey was limited by integrated- 
flux - neither of which are likely. There must, therefore, be 
a distance dependence for column-density sensitivity to av- 
erage galaxies which follow the relationship between Mhi 
and AV. 

Putting in the Mhi sensitivity from Equation 1191 and 
the beam size of Parkes (15 arcmin) gives a column-density 
sensitivity for HIDEEP of: 

Nhi ^ 6.1 x lO^'^dA/pc cm"^ (21) 

Again, this will be affected by the minimum believable 
velocity- width limit. Column-density limits for the three sur- 
veys (HIDEEP, HIPASS BGC and AHISS) are given in Fig- 
ure llSl where it can be seen that HIDEEP has a considerable 
advantage. 

Figure lT^ shows the region of Mhi - Nhi space in which 
each survey is sensitive. HIDEEP extends the region of Mhi 
- Nhi space explored. For instance, neither the HIPASS 
BGC or AHISS would find giant LSB galaxies unless they 
had very high (Mhi / Lb)s or very low velocity widths. 



3.7 Hi properties 

The Hi properties of a sample of HIDEEP sources are given 
in Table 4 (full table available online). Column 1 gives the 
HIDEEP identification for the source. Columns 2 and 3 
give the right ascension and declination from fitting to the 
zeroth-moment map. Columns 4-6 give the noise as mea- 
sured on the part of the spectrum not containing signal (ct) , 
the integrated flux (zeroth-moment, Fhi), and the peak flux 
(•S'peafc), all as measured by mbspect. Columns 7 and 8 give 
the systemic heliocentric velocity (first-moment, Vq) and the 
velocity width at 20% of the peak flux (AV20), measured by 
MBSPECT in the radio- velocity frame (c— ) and converted 
to cz. Column 9 gives the distance in Mpc calculated from 
the CMB rest-frame velocity of the sources, which is approx- 
imately VcMB = V0+278±8 kms"^ for the HIDEEP region 
(the exact correction varies across the field). This does not 
include any correction for bulk-motions (Virgo-centric infaU, 
etc.) nor does it include assigning sources beyond the Cen- 
taurus A group to clusters and assigning a single distance to 
that cluster. Galaxies at less that 1000 kms~^ have been as- 
signed to the Centaurus A group at an assumed distance of 
3.5 Mpc. Column 10 gives the Hi mass calculated from Fhi 
from column 5 and the distance from column 9 using the 
standard equation Mhi = 2.356 x lO^Fif jdlj^^M©, where 
Fhi is in Jy kms^^. 

3.8 Large-scale structure 

The HIDEEP survey region lies in the supergalactic plane 
and there is, therefore, much large scale structure in the sur- 
vey volume. This can be seen in Figure lTTI In particular, the 
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Figure 18. Hi sources and galaxies in NED between 
10,000 kms^i and the HIDEEP bandpass limit of 12,700 kms"!. 
The previously catalogued galaxies with redshifts are shown as 
dots and the HIDEEP sources are shown as filled circles. Abell 
cluster numbers are superimposed on the concentrations of galax- 
ies using positions from Abell, Corwin, & Olowin (1989), these 
have been shifted 0.75 degrees south for the sake of clarity. 



at the north-west of the region. It is possible that the lack 
of correlation is due to the targeting of optical redshift sur- 
veys towards the Abell clusters, thus increasing the number 
of redshifts in those regions out of proportion to their den- 
sity. While Abell 3571 and 1736 are both centred inside the 
HIDEEP volume, at 11,500 kms"^ and 10,500 kms"^ re- 
spectively, Abell 3558 and 3562 are both centred beyond 
14,000 kms"^ (Abell, Corwin, & Olowen 1989; Quintana 
& de Souza 1993). Only the low- velocity tails of these two 
clusters can be seen. 



Figure 17. Large scale structure of galaxies detected in HIDEEP 
(filled circles). These pie-slices show the distribution in R.A. (up- 
per) and Deck (lower), with the angle expanded by a factor of 
three for clarity. Galaxies from NED (crosses) are included for 
comparative purposes, it can be seen the Hi selected sample traces 
approximately the same large scale structure. 



Centaurus A group, the Virgo Southern Extension and the 
Centaurus Cluster can be seen. Near the end of the bandpass 
a number of Abell clusters are found, and more lie beyond 
the outer velocity edge of our survey. 

Figure 1181 shows that there is little correlation be- 
tween the distribution of previously catalogued galaxies with 
redshifts above 10,000 kms~^ and the distribution of the 
HIDEEP sources in this velocity range. Only the Abell 3562 
cluster appears to have a significant number of Hi detec- 
tions and the Hi sources also appear to populate the void 



3.9 Follow-up observations 

Follow up observations of some sources were carried out 
at 21-cm using the Australia Telescope Compact Array 
(ATCA) and optically using the Double Beam Spectrograph 
on the ANU 2.3-m telescope at Siding Spring Observatory. 
The ATCA observations were carried out in the 375-m con- 
figuration in November 1999 and January 2000 and gave 
us high-resolution 21-cm maps of the targeted sources in 
order to accurately determine their positions. Ten out of 
fourteen sources were detected. As the column-density sen- 
sitivity limit for these observations was around 10^^ cm~^, 
substantially higher than the limit for the HIDEEP survey, 
the non-detections do not tell us anything about our survey 
limits. 

Optical spectroscopy has enabled us to positively iden- 
tify ESO 509-G075 as the counterpart of HIDEEP J1335- 
2730 (which was undetected with ATCA) despite a previous 
optical redshift (Quintana et al. 1995) placing it at twice 
the velocity of the Hi source and allowed us to seperate 
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Abell 3558:[MGP94]4312 and 4317 which were too close to- 
gether to be separated by ATCA, identilying 4317 as the 
optical counterpart of HIDEEP J1334-3223. 



4 INFERRED HI COLUMN-DENSITIES 

We have not yet obtained high-resolution Hi images of most 
of our sources so we do not know either their Hi radii or 
column densities. We therefore calculate an inferred column- 
density (N^j) from the Hi mass and the optical (effective) 
radius. Because of our long integration time and low noise, 
we could in principle reach column densities between one 
and two orders of magnitude deeper than previous blind 
surveys (Section I^J. It is of interest to look for evidence, 
however indirect, as to whether this previously unexplored 
parameter space is populated. 

Salpeter and Hoffman (1996) found that thi = (2.34 ± 
0.14) X r_B25 for an optically-selected sample of galaxies. Ob- 
viously r_g25 is not a good measure for the Hi radii of LSB 
galaxies, which may not have a B25 isophote. However, the 
effective radius (re) provides a model and surface-brightness 
independent measure of the optical size of galaxies. In order 
to obtain a relationship between re and thi, it is necessary 
to assume a relationship between rs25 and r^. This obvi- 
ously introduces some unavoidable model-dependency into 
the analysis. As the relationship between vhi and rs25 is 
defined for HSB galaxies, we should define the relationship 
between Ve and r_B25 for a similar sample. There is a fur- 
ther assumption here that the proportionality between 
and rHi found for these HSB galaxies will remain constant 
as we go to lower surface-brightnesses. This may not be the 
case, and would introduce systematic errors into our analy- 
sis, however is certainly a better choice that rs25 for look- 
ing at LSB galaxies. We use the relationship for disc galaxies 
in the ESO-LV catalogue (Lauberts & Valentijn 1989) that 



rs25 = (2.15 ±0.67) x r, 
rni = (5.03 ± 1.59) x rf 



This then gives us: 



(22) 



From this we can calculate the inferred mean column-density 
N?,T as: 



10^ 



(23) 



where Rhi is the radius in pc, calculated from vhi using 
the distance to the source. 

Figure IT^ shows these inferred column-densities plotted 
against their velocity widths. Above the survey limit (dashed 
line), the volume sampled depends only on Hi mass (or, 
more precisely, on peak flux) and not on column-density, 
as shown earlier in Figure If low column-density galaxies 
were common, we would expect to see galaxies at all column- 
densities in this plot. However, this is clearly not the case 
- there is a large gap between the lowest column-density 
galaxies found and our sensitivity limit. 

That this gap is real and not an artefact of our method is 
shown in Figure 1^ Here the sample has been split into four 
sections, according to the beam-filling factor {AQ/Qi,eam) of 
the sources, and the column-density limits calculated with 
this included. The top left panel shows that we are detecting 
galaxies up to the column density limit of the survey for this 
range in beam filhng factors. These are galaxies with small 
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Figure 19. Inferred column-densities of HIDEEP sources against 
velocity widths. The limit shown is for the 18 mjy completeness 
limit for HIDEEP and sources filling the beam. 
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Figure 20. Inferred column-densities of HIDEEP sources with 
diff'erent beam filling factors against velocity widths. The limits 
shown are for the 18 mJy HIDEEP completeness limit, the dotted 
line indicating the limit for the lower filling factor and the solid 
line the limit for the higher filling factor in each subgraph. 



beam-filling factors, and we can see that the inferred column 
density limits for these galaxies are quite high (cf. Eqn B6) . 
The top right panel shows the limits for galaxies with slightly 
larger beam filling factors. Galaxies are detected with col- 
umn densities close to the survey limit for this range in beam 
filling factors, but we can see that they do not approach the 
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Figure 21. Distribution of inferred column-densities of HIDEEP Figure 22. Relationship between effective radius and Hi 

sources. Tlie dashed line indicates a Gaussian with a mean of flux. The dashed line shows the expected slope for N'^^ = 

20.65 and a scatter of 0.38. It can be seen that this is a fairly ]^q20.65 2^ 
good description of the distribution. 



limit as closely as in the previous panel. This is an indication 
that in this sample galaxies with larger beam filling factors 
(and hence lower column densities; cf. Sect. 1) are rarer. This 
is confirmed by the two bottom panels. These show the lim- 
its for galaxies that come close to filling the beam (and thus 
have the lowest column densities). We see that despite the 
low limits, there are only two galaxies detected. This thus 
indicates an absence of low-column density galaxies, despite 
their potential detectability. 

Using the numbers of galaxies found within the over- 
lap of the full-sensitivity Hi survey and the optical survey, 
we calculate that, to 95% confidence, low column-density 
galaxies make up less than 21% of galaxies with Hi masses 
between 10* and ICFMq (13 galaxies), less than 6% between 
10^ and 1O^°M0 (52 galaxies), and less than 19% between 
10^° and 10"Mq (14 galaxies). 

If our estimate of Rhi were out by more than a factor 
of three, we would expect to see a lower limit to our col- 
umn densities parallel to the plotted completeness limit. As 
the lower limit appears flat with respect to velocity width, 
it appears that there are indeed no large very low column- 
density galaxies in our sample. If such objects are present 
in the local Universe they must, therefore, be rare. We dis- 
cuss possible reasons for this absence elsewhere (Disney & 
Minchin 2003). 

The distribution of column densities appears to be 
the same at all surface brightnesses. A two-dimensional 
Kolmogorov-Smirnov test (Peacock 1983) shows that the ob- 
served distribution is indeed consistent with a distribution 
having the same column-density distribution as HIDEEP 
(Fig. at every surface brightness. The distribution 

shown in Figure 1211 is well described by a Gaussian with 
log(iVijj/cm"^) = 20.65 ± 0.38. Figure l22l shows the mea- 
sured data which gives rise to a constant column-density: 



a relationship between the effective radius and the Hi flux. 
The dashed line on this graph is for a constant N'^j of 10^°'^^ 
cm~^, and it can be seen that this matches the data well. 
This apparent constancy of column density is unexpected, 
we would expect to see a fall-off towards lower surface- 
brightness, as observed by de Blok et al. (1996) for a sample 
of LSB galaxies observed with the Westerbork Synthesis Ra- 
dio Telescope (WSRT). 

Swaters et al. (2002) observed a number of optically- 
selected galaxies with the WSRT, showing that the relation- 
ship between Hi mass and diameter was much tighter when 
true Hi diameters were used than when they were inferred 
from the optical radii. It is likely that much of the scatter 
in our relationship has been similarly introduced by our use 
of optical radii. Observations of the true Hi diameters may 
well show a similar trend to that seen by de Blok et al. 

In order to compare our sample with the literature, we 
have taken the LSB galaxy catalogue of Impey et al. (1996; 
ISIB96). This survey provides effective surface-brightnesses 
and radii for all the galaxies found, across a wide range 
of surface brightness, and provides Hi mass measurements 
for a subsample of 190 galaxies. We therefore calculate the 
column density from this data in exactly the same manner 
as for HIDEEP. The Impey et al. survey was carried out in 
_B-band, we therefore convert /if to ^f using the average 
colour for disc galaxies from de Jong 1996, B — R — 1.1. 

Figure 1231 shows the surface-brightness - column- 
density distribution for both surveys. HIDEEP galaxies are 
indicated by solid triangles, ISIB96 galaxies by open circles. 
It can be seen that there is a clear trend in the ISIB96 data, 
which is similar to that seen by de Blok et al. (1996). The 
surface-brightness distribution of the ISIB96 galaxies is dif- 
ferent to that of HIDEEP, due to the method in which the 
galaxies have been selected. In order to compare the sur- 
veys we have therefore re-sampled the ISIB96 data to give 
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Figure 23. Comparison between tlie surface-brightness — 
column-density distribution for HIDEEP (triangles) and Impey 
et al. (1996), converted to il-band assuming B — R = 1.1 (open 
circles). Error bars have been omitted for the sake of clarity. 



the same surface-brightness distribution (in 0.5 mag bins) 
as HIDEEP. This was carried out ten times, and the sub- 
samples compared with HIDEEP using the two-dimensional 
Kolmogorov-Smirnov test. Half of the subsamples were not 
significantly diff'crcnt at the 95% level. 

There are a number of systematics that could affect this 
comparison. Our measure of re could be systematically too 
small towards lower surface-brightnesses: this would lead to 
both Nffj and /le being higher and so would act to destroy 
such a relationship between Nhi and surface-brightness. 
The colours will not be exactly the same for all the galax- 
ies, introducing errors into the conversion of the ISIB96 
relationship to 7J-band. The low column-density galaxies 
i^Hi < 10^" cm~^) in ISIB96 are generally large galaxies, 
with fcff > 20". Using our scaling between thi and r^ff, 
this gives them Hi sizes larger than the Arecibo beam. It is 
therefore likely that some of the hydrogen in these galaxies 
was missed in the pointed observations, leading to lower in- 
ferred column-densities than is truly the case. We cannot, 
therefore, say that the HIDEEP data is inconsistent with 
ISIB96. If there is a relationship between surface-brightness 
and Nhi in our data, however, it is very weak indeed. 



5 DISCUSSION, CONCLUSIONS AND 
FUTURE WORK 

The existence of a large population of LSB galaxies could 
have a significant impact on several areas of astronomy, in 
particular on measurements of the luminosity and mass den- 
sity of galaxies and on theories of galaxy formation and evo- 
lution. However, LSB objects are difficult to detect in the op- 
tical since, by definition, they are hard to discriminate from 
the background sky. Hence, alternative techniques need to 



be considered for determining whether there exists a signifi- 
cant population of extragalactic objects with optically LSB. 
Searches for galaxies using the 21-cm line of atomic hydro- 
gen provide one such alternative method. If LSB galaxies 
have similar hydrogen content to galaxies of higher surface 
brightness, then we might expect that a 21-cm survey will 
not be biased against LSB galaxies. 

However, the detection of a galaxy in a 21-cm survey is 
a function of both its Hi column density and its Hi mass. 
If a galaxy has a column density below the column density 
limit of the survey then it will never be detected whatever 
its total Hi mass. Assuming LSB galaxies have similar Hi 
content to brighter galaxies but that this is distributed over 
a larger surface area, then we would expect LSB objects to 
have lower Hi column densities. Previous Hi surveys have 
not generally had long enough integration times to search 
to sufficiently low column densities to draw definite conclu- 
sions about the cosmic prevalence of gas-rich LSB galax- 
ies. HIDEEP with its 9000s beam" ^ integration time is the 
first survey with sufficient column density sensitivity to be 
able to place interesting limits on the existence and size of 
a low column-density/LSB population. Two major conclu- 
sions from HIDEEP have been presented here. 

Firstly, all of the sources found in HIDEEP appear to 
be associated with am optical counterpart on our deep UK 
Schmidt R-band data. In other words, we have not found 
an intergalactic Hi cloud down to our observational column 
density limit of N/jz/AV = 2.1x10^*' cm"^ which corre- 
sponds to an Nhi of 4.2x10^^® cm~^ for a typical galaxy 
velocity-distribution (AV=200 kms"^) and 3.2xl0" cm"^ 
for a typical QSOAL dispersion (15 kms~^). Wherever neu- 
tral hydrogen is found it is accompanied by a visible popu- 
lation of stars above our surface-brightness limit of 26.5 R 
magnitudes per square arcsecond (~27.5 in B). 

Secondly, if we infer the Hi sizes of the HIDEEP detec- 
tions from their optical sizes, then we can derive Hi column 
densities for them. These derived column densities are all 
>10^'' cm~^, more than an order of magnitude above our 
observational limit. Assuming that our method of inferring 
Hi size from optical size is robust, then this result implies 
that there is no significant population of galaxies with Hi 
column densities <10^'^ cm~^. We are currently undertak- 
ing a VLA and ATCA survey of all the HIDEEP sources in 
order to measure their column densities directly. This will 
enable us to test this result without having to infer Hi size 
from optical size. 

Possible physical explanations for this intriguing second 
result are discussed in detail elsewhere (Disney & Minchin 
2003). lonisation is unlikely to be responsible since the inter- 
galactic radiation field locally is at least an order to magni- 
tude too small (Scott et al. 2002). Another possibility, raised 
in the early days of 21-cm astronomy, is 'freezing out' - i.e. 
that the 21-cm spin temperature of low column density ob- 
jects may fall to the cosmic background temperature, ren- 
dering such clouds invisible in emission. 

The cosmic significance of gas-rich low-surface bright- 
ness galaxies will be the subject of a subsequent paper based 
on this same data (Minchin et al. 2003, in preparation). 
However, the fact that we have not found a single galaxy in 
Hi which cannot be seen in our optical data suggests that 
there cannot be large numbers of gas-rich extremely LSB 
galaxies or intergalactic clouds. 
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If our second result is confirmed by the follow-up VLA 
and ATCA observations, then this will imply that Hi sur- 
veys need not have column density limits much lower than 
^10^° cm~^, since few galaxies would appear to have col- 
umn densities lower than this. In particular, this would im- 
ply that previous large-area, shallow surveys, e.g. HIPASS 
and HIJASS, will not have missed many low-column density 
objects. However, long integrations using telescopes with 
smaller beams would improve the statistics on the numbers 
of low mass, low column-density sources. 
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APPENDIX A: DERIVATION OF THE 
SURFACE-BRIGHTNESS - COLUMN-DENSITY 
RELATION 

For any given area the Hi surface density and the optical 
surface brightness arc related by: 

^Hi (Mopc-^) = X Es (LqPC-^) (Al) 

where T,hi is the Hi surface density and S_b is the B- 
band optical surface brightness, averaged over the region, 
and Mhi and Lb are the Hi mass and the i?-band lumi- 
nosity within the same region. As 1 Mq pc~^ is approxi- 
mately equal to an Hi column density of 10^"'^ cm~^ and 1 
Lq pc^^ is approximately equal to a surface-brightness of 
27.05 _Bmagarcsec~^, this gives the scaling relationship; 

Nhi ^ 10^°-^ (^) lo'"-*'^^-''— » (A2) 

where Nhi is the Hi column density in units of cm^^ and 
A*mean is the average optical surface brightness in units of 
magarcsec"^ taken over the same area as Nhi- This can be 
re-written as: 
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/in 



log {Nh 



(A3) 



This relation can be adapted to relate the central 
surface-brightness of a galaxy to its average Hi column den- 
sity if certain assumptions are made about the size of the 
Hi disc. Cayatte et al. (1994) found that Rhi ~ 1.7i?25- 
This scaling will obviously not hold for LSB galaxies, some 
of which may not even have a /is =25 isophote, but for a 
Freeman's law galaxy _R25 =3.1 scale lengths and it there- 
fore seems reasonable to assume that Rhi = 3.1 x 1.7 = 5.25 
scale lengths. Assuming also that Mhi/Lb = O.3M0/L0 
(average value from Roberts & Haynes 1994), we get: 



^0 ^ 2.5 (28.95 - log [Nhi)) 



(A4) 



which can be used to work out an approximate equivalent 
central surface-brightness limit for Hi surveys. 



APPENDIX B: COLUMN-DENSITY 
SENSITIVITY 

The signal entering the receiver is measured in terms of the 
antenna temperature Ta where 



kTA = S,,D' 



(Bl) 



Sv is the strength of the source in flux units and D the 
dish-diameter. For a signiflcant detection Ta should exceed 
the uncertainty in the system power by some signal-to-noise- 
ratio cr, i.e. 



(B2) 



which is the usual 'antenna equation' derived from the band- 
width theorem. 

In the Rayleigh- Jeans regime, surface brightness is con- 
ventionally expressed in terms of the brightness temperature 



Tb = — 



2k An 



(B3) 



where Afl is the solid angle of the source. For Hi galaxies 
quantum mechanics yields 



^ = 1.8 X w'^Tb 



(B4) 



where Nhi is the column-density in cm~^ and AV is the 
velocity width of the line. 

The antenna equation IIB2l l can now be rewritten as 



^> 1.8x10- ' 



ciTs, 



2D^AQ, ^tobsAu,^ 



(B5) 



where we have replaced Ai/ with Ai/ch because we are here 
interested in a peak-flux limited survey. 

If we now substitute in the beam size, Qb = 1.13(A/Z))^, 
and the channel velocity-width AVch = (c/urest) x Avch, we 
get 



Nhi 



> 5.5 X 10 



16 ^'fc 



AV ■ An VtobsAV^H 

So, if the source fills the beam, i.e. ilb/AQ = 1 



Nhi 
AV 



> 5.5 X 10" 



; ^ Tsys 



(B6) 



(B7) 



which agrees with Equation |S| and is indeed independent of 
telescope diameter. For HIDEEP, T^ys = 26K, Ubs = 9000s, 
AVch = 13.2 kms~^ and if ct = 5 

Nhi 



> 2.1 X 10' cm~7kms~ 



(B8) 



AV 

For an optimally smoothed survey limited only by the 
signal-to-noise of the total flux of the galaxy rather than 
the flux in a single channel, one would replace AVch by 
(AV/c)v2i in Eqnl^and reach (Disney & Banks 1997) 



'^"'^ >1.8xl0-ar„ 



(Nhi\ 
\AV ) 



AVtcb 



(B9) 



APPENDIX C: CALCULATING MINIMUM 
COLUMN-DENSITY SENSITIVITY FROM THE 
LOWEST FLUX SOURCES ACTUALLY 
DETECTED 

At any distance a source that is just detectable because of 
its Hi mass, and which just flUs the beam, will have the 
minimum column-density the survey is capable of detecting: 

20 Mhi 



Nhi = 5 X 10^ 



-KO^dlc 



(CI) 



where 9 is angular diameter of the source in radians, and 
dpc its distance in parsecs. By a well known relation Mhi is 
related to flux Fhi by: 



Mhi = 2.356 x lO^F^/d^ 



Mpc 



(C2) 



and therefore we can substitute this into the equation above 
to give 



Nhi 



2.356 X lO^fff/difpc 



■KO^dlc 



therefore 



iVJ?r = 4.5 X 10'° 



i^HI 



(C3) 



(C4) 



where 6' is the source diameter in arcminutes. Putting 9' 
equal to the beam diameter, and Fhi to the lowest flux 
measured in the survey will yield the column-density limit 
of the survey at the velocity width of that source. Of more 
interest is putting in the lowest value of Fhi/AV, which 
allows the column-density limit to be found at all velocity 
widths: 



(NhiV 
\AV J 



.4.5x,0»(^) 



AV 

which is Ean ll2l in the text 



mm 2^ 
fl72 
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